A microchip optomechanical accelerometer 
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The monitoring of accelerations is essential for a variety of applications ranging from inertial navigation to 
consumer electronics [Tj. The basic operation principle of an accelerometer is to measure the displacement 
of a flexibly mounted test mass; sensitive displacement measurement can be realized using capacitive [2 3], 
piezo-electric |4|, tunnel-current |5|, or optical |6-9| methods. While optical readout provides superior dis- 
placement resolution and resilience to electromagnetic interference, current optical accelerometers either do 
not allow for chip-scale integration |6| or require bulky test masses |7 8|. Here we demonstrate an optome- 
chanical accelerometer that employs ultra-sensitive all-optical displacement read-out using a planar photonic 
crystal cavity [10] monolithically integrated with a nano-tethered test mass of high mechanical g-factor 1 1 1 1. 
This device architecture allows for full on-chip integration and achieves a broadband acceleration resolution of 
10 flg/\/Hz, a bandwidth greater than 20 kHz, and a dynamic range of 50 dB with sub-milliwatt optical power 
requirements. Moreover, the nano-gram test masses used here allow for optomechanical back-action [ 12 1 in the 
form of cooling 1131 or the optical spring effect I14II15I . setting the stage for a new class of motional sensors. 
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Due to the rapid development of silicon micro machining 
technology, MEMS accelerometers have become exceedingly 
popular over the last two decades [1]. Evolving from airbag 
deployment sensors in automobiles to tilt-sensors in cameras 
and consumer electronics products, they can now be found 
in a large variety of technological applications with very di- 
verse requirements of their performance metrics. While sen- 
sors for inertial navigation systems require low noise levels 
and superior bias stability [ 16 1, large bandwidth is crucial for 
sensors in acoustics and vibrometry applications. However, 
there is a fundamental tradeoff between noise performance 
and bandwidth which can be understood from the basic op- 
eration principle of an accelerometer, illustrated in Fig. [Tk. 
When subjected to an acceleration a((o) at frequency CO, a 
mechanically compliant test mass experiences a displacement 
x(co) = x((o)a(a>) proportional to the mechanical susceptibil- 



ity X ((0) = cof n -co l + i°^-. Here, (O m = 2nf m - 
the (angular) resonance frequency of the oscillator and Q„ 
its mechanical Q-f&c\.ox (see the plot of |#(fi>)| in Fig.fTh for 
Q,„ = 10). Usually, accelerometers are operated below their 
fundamental resonance frequency co m , where %{co) « l/ft>, 2 
exhibits an almost flat frequency-response. This naturally 
leads to a tradeoff between resolution and bandwidth, since 
the large resonance frequency required for high-speed opera- 
tion results in vanishingly small displacements. As a result, 
the performance of the displacement sensor constitutes a cen- 
tral figure of merit of an accelerometer. 

In a cavity optomechanical system, a mechanically compli- 
ant electromagnetic cavity is used to resonantly-enhance read 
out of mechanical motion [17] (canonically, the motion of the 
end mirror of a Fabry-Perot cavity). Such systems have en- 
abled motion detection measurements with an imprecision at 
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or below the standard quantum limit (SQL) [18-20], corre- 
sponding to the position uncertainty in the quantum ground- 
state of the mechanical object. Clever quantum back-action 
evading techniques [21 1 aside, only for an ideal cavity sys- 
tem (no parasitic losses) can the actual displacement sensitiv- 
ity reach the SQL due to fluctuating radiation pressure forces 
arising from shot noise of the probe light [22|. The aver- 
age radiation pressure force, on the otherhand, can be quite 
large in micro- and nano-scale optomechanical devices, and 
offers the unique capability to control the sensor bandwidth 
via the optical spring effect lfl4l [T31 and the sensor's effec- 
tive temperature via passive damping |[T2l or feedback cold- 
damping mum. 

In this work, we utilize an integrated silicon-nitride (SiN) 
zipper photonic crystal optomechanical cavity 1 10 ] to pro- 
vide shot-noise-limited read out of mechanical motion with 
imprecision at the SQL, enabling high-bandwidth and high- 
resolution acceleration sensing. The resolution of an ac- 
celerometer can be quantified by a noise-equivalent acceler- 
ation, NEA = a /a 2 ,, + oL + a 2 dd in units of g/v 7 ^ (1 g = 



h ^ "det 

9.81 m/s 2 ). The first term in the NEA is due to thermal Brow- 



nian motion of the test mass (see appendix III [24| and is 
given by, 



0th 



l4k B Tca m 



mQ„ 



(1) 



while the remaining two terms arise from the aforementioned 
displacement readout noise («det) and added noise (back- 
action) onto the test mass due to the act of measurement (a a dd> 
see appendix [T~4| >. Fundamental to minimizing the NEA is a 
reduction in the intrinsic thermal noise, a X h, which according 
to equation ([TJi, requires one to maximize the mass-Q prod- 
uct at a given co m . In most commercial accelerometers, the 
g-factor is relatively low, which demands large test masses 
for high resolution. In contrast, in the zipper cavity devices 
presented here, we use nano-tether suspension of a nano- 
gram test mass to yield high intrinsic mechanical Q-factors 




FIG. 1: Overview of the accelerometer design, a, Canonical example of an accelerometer. When the device (blue frame) experiences a 
constant acceleration a, a test mass m undergoes a displacement of x = ma/k. b, Frequency response \% (fi))| of an accelerometer on a log-log 
plot featuring a resonance at f m = sjk/m/2iz with Q m = 10. c, False-colored SEM-image of a typical optomechanical accelerometer. A test 
mass of size 150 flm x 60 jUm x 400 nm (green) is suspended on highly stressed 150 nm wide and 560 flm long SiN nano-tethers, which allow 
for high oscillator frequencies (> 27 kHz) and high mechanical g-factors (> 10 ). On the upper edge of the test mass, we implement a zipper 
photonic crystal nanocavity (pink). The cross-shaped cuts on the test mass facilitate undercutting the device, d, Zoom-in of the optical cavity 
region showing the magnitude of the electric field |E(r)| for the fundamental bonded mode of the zipper cavity. The top beam is mechanically 
anchored to the bulk SiN and the bottom beam is attached to the test mass, e, Schematic displacement profile (not to scale) of the fundamental 
in-plane mechanical mode used for acceleration sensing, f, SEM-image of an array of devices with different test mass sizes. 



(1 —2 x 10 6 ), and strong thermo-optomechanical back-action 
to damp and cool the thermal motion of the test mass. 

Figure [Tt shows a scanning-electron microscope image 
of the device studied here, with the test mass structure and 
nano-tethers highlighted in green. The fundamental in-plane 
mechanical mode of this structure is depicted in Fig. [it 
and is measured to have a frequency of f m — 27.5 kHz, in 
good agreement with finite-element-method simulations from 
which we also extract a motional mass of m — 10 x 10~ 12 kg. 
The measured mechanical Q-f&c\ar is Q,„ = 1.4 x 10 6 in vac- 
uum (see appendix [Gb, which results in an estimated a t h = 
1.4 /xg/^/Hz. The region highlighted in pink corresponds to 
the zipper optical cavity used for monitoring test mass mo- 
tion, a zoom-in of which can be seen in Figure [TJl. The cav- 
ity consists of two patterned photonic crystal nanobeams, one 
attached to the test mass (bottom) and one anchored to the 
bulk (top). The device in Fig. [it is designed to operate in 
the telecom band, with a measured optical mode resonance at 
X = 1537 nm and an optical g-factor f Q = 9,500. With 
the optical cavity field being largely confined to the slot be- 
tween the nanobeams, the optical resonance frequency is sen- 
sitively coupled to relative motion of the nanobeams in the 
plane of the device (the i-direction in Fig. [TJ;). A displace- 
ment of the test mass caused by an in-plane acceleration of the 
supporting microchip can then be read-out optically using the 
setup shown in Fig. [2k, where the optical transmission through 
the photonic crystal cavity is monitored via an evanescently- 
coupled fiber taper waveguide [25 1 anchored to the rigid side 
of the cavity. Utilizing a narrow bandwidth (< 300 kHz) laser 
source, with laser frequency detuned to the red side of the cav- 
ity resonance, fluctuations of the resonance frequency due to 
motion of the test mass are translated linearly into amplitude- 
fluctuations of the transmitted laser light field (see inset in 



Fig. [2k and appendixE). A balanced detection scheme allows 
for efficient rejection of laser amplitude noise, yielding shot- 
noise limited detection for frequencies above ~ 1 kHz. 

Figure Wr> shows the electronic power spectral density 
(PSD) of the optically transduced signal obtained from the 
device in Fig. [it. The cavity was driven with an incident 
laser power of P, n = 116 /J.W, yielding an intracavity photon- 
number of fa 430. The two peaks around 27.5 kHz arise from 
thermal Brownian motion of the fundamental in- and out- 
of-plane mechanical eigenmodes of the suspended test mass. 
The transduced signal level of the fundamental in-plane reso- 
nance, the mode used for acceleration sensing, is consistent 
with an optomechanical coupling constant of goM = 2n x 
5.5 GHz/nm, where goM = d(0 /dx is defined as the optical 
cavity frequency shift per unit displacement. The dotted green 
line depicts the theoretical thermal noise background of this 
mode. The series of sharp features between zero frequency 
(DC) and 15 kHz are due to mechanical resonances of the an- 
chored fiber-taper. The noise background level of Fig. [2J5 is 
dominated by photon shot-noise, an estimate of which is indi- 
cated by the red dotted line. The cyan dotted line in Fig. [2b 
corresponds to the electronic photodetector noise, and the pur- 
ple dashed line represents the sum of all noise terms. The 
broad noise at lower frequencies arises from fiber taper mo- 
tion and acoustic pick-up from the environment. The right- 
hand axis in Fig. [2b quantifies the optically transduced PSD 
in units of an equivalent transduced displacement amplitude 
of the fundamental in-plane mode of the test mass, showing a 
measured shot-noise-dominated displacement imprecision of 
4 ftn/vHz (the estimated on-resonance quantum-back-action 
displacement noise is 23 fm/VHz, and the corresponding on- 
resonance SQL is 2.8 fm/\/Hz; see appendix '. 
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At this optical power the observed linewidth of the mechan- 
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FIG. 2: Experimental system and noise data, a, Laser light used 
to probe the zipper cavity motion is split with a beamsplitter; the 
signal arm is sent through a fiber polarization controller (FPC) and 
a fiber taper, which is coupled to the optical cavity, while the other 
arm is sent directly to a balanced photo-detector (BPD). Variable op- 
tical attenuators (VOA) in each arm balance the powers, and a power 
meter (PM) is used to calibrate the probe power. The BPD signal is 
sent to a proportional-integral controller (PI) - locking the laser half 
a linewidth red-detuned from the optical resonance. The sample is 
mounted on a shake table comprised of a shear piezo. Transduced 
accelerations are measured using either an electronic spectrum ana- 
lyzer (ESA) or a lock-in amplifier, b, The left axes show an optical 
power spectral density (PSD) plot of the BPD signal showing me- 
chanical modes at 27.5 kHz (green). The right axis shows the equiv- 
alent displacement noise. The tone at 26 kHz (orange) is transduction 
of a tone applied to the shear piezo corresponding to an acceleration 
of 38.9 mg. The dashed and dotted lines are theoretical noise levels 
for shot noise (red), detector noise (cyan), thermal noise (green), and 
the total of all noise contributions (purple). The inset is a time trace 
of the transduction of an applied acceleration of 35.6 mg at 25 kHz. 



ical mode is w 2 Hz, roughly 100 times larger than the low 
power linewidth. As modeled in appendix [H] the measured 
mechanical damping is a result of radiation pressure dynam- 
ical back-action, enhanced by slow thermo-optical tuning of 
the cavity which provides the necessary phase-lag for efficient 
velocity damping. Damping of the mechanical resonance is 
typically used to reduce the ringing transient response of the 
sensor when subjected to a shock input [26|. In contrast to 
conventional gas-damping employed in MEMS sensors ll27l . 
optomechanical back-action damping also cools the mechani- 
cal resonator [ 13 1. The measured effective temperature of the 
fundamental in-plane mode of the test mass, as determined 



from the area under the 27 kHz resonance line in Fig. [2b, is 
r e ff f=a 3 K. This combination of damping and cooling keeps 
the ratio of T e ff/Q m fixed, and does not degrade the thermally- 
limited acceleration resolution of the sensor. 

In order to carefully calibrate the accelerometric perfor- 
mance of the device, the sample is mounted onto a shake table 
driven by a shear piezo actuator (see appendixfFli. Applying a 
sinusoidal voltage to the piezo results in a harmonic accelera- 
tion fl(o), and thereby a modulation of the transmitted optical 
power. The optical power in the modulation sideband is given 
by (see appendix [El 



P m {co) = (1 - T d ) ^Pingou \x(co)a(co)\, 

(On 



(2) 



where Q is the optical (Mactor (= 9,500), (0 is the optical 
resonance frequency, T c / is the relative cavity transmission on 
resonance (= 0.88), and the laser is half a linewidth detuned. 
The narrow tone at 26 kHz in Fig. [2j3 (orange) arises from an 
applied rms-acceleration of a lms = 38.9 mg, calibrated using 
two commercial accelerometers mounted on the shake table 
(see appendix [Fl, From the signal-to-noise-ratio of this cali- 
bration tone we estimate a± = 2.0 jJ-g/y/Hz, comparable to 
the theoretical value of a^ = 1.4 /xg/\/Hz. For a driving tone 
at 10 kHz, we measure a m ; n s=s 10 /ig/-\/Hz, limited in this 
case by photon shot noise. The dynamic range over which the 
sensor is linear at a drive frequency of 10 kHz has also been 
measured (see appendix 111), and is found to be > 49 dB (up to 
~ 10 g accelerations, limited by the maximum output voltage 
of the piezo shaker drive electronics). 

Figure [3b shows the demodulated photodiode signal nor- 
malized to the applied acceleration as a function of drive fre- 
quency, corresponding to the frequency dependent accelera- 
tion sensitivity of the zipper cavity (the inset of Fig. [3k shows 
data from the commercial accelerometers used to calibrate 
the applied acceleration). The dashed red line is the theoret- 
ical calculation of the sensitivity without fit parameters and 
shows excellent agreement. The sharp Fano-shaped features 
for lower frequencies can again be attributed to mechanical 
resonances of the fiber-taper waveguide. The broad region of 
apparent higher-sensitivity around 15 kHz is due to an under- 
estimate of the applied acceleration arising from an acoustic 
resonance of the shake table. 

The calibrated frequency-dependent NEA, shown in 
Fig. [3b, is obtained by normalizing the ESA noise spectrum 
(Fig! T5V)) by the sensitivity curve (Fig. [3k). Between 25- 
30 kHz the resolution is limited by the thermal noise of the os- 
cillator, while from 5-25 kHz shot-noise limits the resolution 
to w 10 /ig/\/Hz. For frequencies lower than 5 kHz, motion 
of the fiber-taper waveguide and the environment contribute 
extra noise. The sharp Fano-shaped feature at 27 kHz arises 
from interference with the fundamental out-of -plane mode of 
the test mass. The dashed red curve corresponds to a theoret- 
ical estimate of the NEA which shows good agreement. The 
dashed green line is the fundamental thermal sensing limit. 

The device platform demonstrated here straightforwardly 
allows for further reduction of the NEA. For instance, a t h can 
be reduced further by increasing the test mass m. In a prelim- 
inary study, we have fabricated a series of devices with test 
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FIG. 3: Frequency-dependence of sensitivity and resolution, a, Sensitivity curve as function of frequency, obtained by driving the shear 
piezo with a sinusoidal voltage and measuring the amplitude of the resulting voltage modulation of the BPD signal using a lock-in amplifier. 
The dashed red line corresponds to the theoretical expectation for the sensitivity without fit parameters. The inset shows data from commercial 
accelerometers also attached to the shake table (blue and green curves), which are used for calibrating the applied acceleration, b, Frequency- 
dependent noise-equivalent acceleration (NEA) of the device quantifying its broadband-resolution. The plot is obtained by taking the PSD in 
Fig.pJ? and normalising it by the sensitivity curve in a. The dashed red line depicts the theoretical expectation for the NEA given shot-noise 
and thermal noise limitations. The green dashed curve corresponds to the thermal noise (a t h)- 



masses ranging from 100 x 10~ 15 kg to 35 x 10~ 12 kg and 
recorded their mechanical frequency and Q-f&clor:. Figure Ht 
depicts the calculated a t h versus the mechanical frequency of 

3/2 

the studied devices, which roughly scales with a t h x fiW 
(green line). Adding mass alone also results in a reduction 
of the sensor bandwidth; however, by scaling the number 
of nano-tether suspensions with the test mass size (see Fig- 
ures Rb and c) the bandwidth can be kept constant. More- 
over, as shown in the inset of Fig. Hk, we have found that 
adding nano-tethers does not result in a degradation of the 
mechanical 2-factor. Simultaneously scaling the width of the 
test mass and the number of nano-tethers by a factor of 100 
from the device shown in Fig. lie, to a mass of m = 10~ 9 kg, 
should reduce the thermal NEA to ~ 150 ng/v^Hz while 
maintaining a sensor bandwidth of 25 kHz. Critically, for 
gOM = 2nx 100 GHz/nm as measured in previous zipper cav- 
ity structures [ 10], the optical input power required to reach 
this resolution across the entire sensor bandwidth is still sub- 
milliwatt (~ 850 jiW). 

With a demonstrated acceleration resolution on the order of 
a few /xg/\/Hz and a bandwidth above 25 kHz, the zipper cav- 
ity device presented here shows performance metrics orders 
of magnitude better than other optical accelerometers J6] |7) 
and comparable to the best commercial sensors [28 1. These 
devices, formed from a silicon chip, also allow for the inte- 



gration of electrostatic tuning -capacitors ||29ll , fiber-coupled 
on-chip waveguides [8 1, and on-chip electronics, all of which 
enables convenient, small form-factor packaging, and elim- 
inates the need for expensive tunable lasers. In addition, 
nanoscale optomechanical cavities such as the zipper cavity 
studied here, offer the unique resource of strong radiation- 
pressure back-action. The optical spring effect, for example, 
allows for dynamic tuning of the mechanical resonance fre- 
quency, which can increase the low-frequency displacement 
response (inverse quadratically with frequency) and decrease 
thermal noise (with the square root of frequency). Similar 
zipper cavity devices have shown low power (sub-mW) op- 
tical tuning of the mechanical resonance frequency over 10's 
of MHz (> 200% of co m ) into a regime where the mechanical 
structure is almost entirely suspended by the optical field [ 10 1. 
Also, as demonstrated here, back-action cooling provides a re- 
source to damp the response of the oscillator without compro- 
mising the resolution. Combining all of these attributes should 
allow not only for a new class of chip-scale accelerometers, 
but other precision displacement-based sensors of, for exam- 
ple, mass, force, and rotation. 
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Appendix 
Appendix A: Oscillator susceptibility 



The oscillator susceptibility #(©) given above follows from the differential equation of the harmonic oscillator: 
Transforming to Fourier space, this reads 



2 

mx + mjx + m(O m x = F app i . ( A 1 ) 



,2„ i ,•„,„„ i ™2 



^appl(<» 



cox + icoyx + a>f n x = — — . (A2) 



With f ap pi(o)/m = a ap pi, this yields the accelerometer response 

x(co) = lOKppiO) 



1 



^aappi(w). (A3) 



This function has the following properties: 



X(co m ) = -i-§ = -iQmX{0), (A5) 

X(o)>o) m ) <* — _. (A6) 

For the device studied here with ©„, = 27.5 kHz, this gives an acceleration sensitivity of j(0) = 329 pm/g with g = 9.81 m/s 2 . 

Appendix B: Sample fabrication and design 

The presented accelerometer structures are defined in a 400 nm thick silicon nitride (SiN) layer formed on top of a 500 jUm 
thick single-crystal silicon wafer. The SiN is stiochiometric and is grown in LPCVD under conditions that allow for large internal 
tensile stress (a = 800 MPa). The accelerometer structures comprising the test mass, the support nano-tethers, and the zipper 
cavity are defined in a single electron-beam lithography step. The mask is transferred into the SiN layer using ICP/RIE dry- 
etching in a SF6/C4F8 plasma. Resist residues are removed in a combination of heated Microposit 1 165 remover and Piranha 
solution (3:1 H2SO4 : H2O2) at 120° C. The structures are undercut by anisotropic wet-etching in 70° C hot KOH and cleaned 
in a second Piranha etching step. Critical point dying in CO2 avoids collapsing of the zipper cavities. 

The optical and mechanical structures are designed using finite-elements simulations performed in COMSOL Multiphysics 
(http://www.comsol.com/). 

Appendix C: Optical spectroscopy 

The sample is optically coupled via a near-field probe consisting of a tapered optical fiber. The tapered fiber is brought in 
optical contact with the device using attocube nanopositioners. Aligned in parallel to the zipper nano-beams, the fiber taper is 
mechanically anchored on the struts attached to the rigid side of the zipper cavity. Launching light from a NewFocus Velocity 
tunable external-cavity diode laser into the fiber taper and monitoring the taper transmission then allows us to do resonant 
coherent spectroscopy of the cavity mode. Technical amplitude noise of the laser (~ 10 dB above the shot-noise level) is 
suppressed by a balanced detection scheme using a Newport 21 17 balanced photodetector that features ~ 20 dB common-mode 
noise rejection. 

Appendix D: Transmission function of side-coupled open cavity 

In order to calculate the intensity transmission profile T(co) of a photonic-crystal resonator side-coupled by a fiber-taper 
waveguide, we start from the equation of motion of a, the annihilation operator of the cavity field: 

dfl / K\ „ /fQ„ [Kg ^ /A „ N 

— = -[iA+-)a + J — dm + VKiOi + J — a-. (A7) 
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FIG. Al: Example transmission curve of the zipper cavity. The curve is obtained by scanning an external cavity diode laser across the 
cavity resonance at Ay = 1537.36 nm while monitoring the fiber taper transmission. The dip exhibits an optical g-factor of Q„ = 9,500 and a 
transmission dip on resonance of Tj = 0.88. 



Here, A = 0| — <a c is the laser-cavity detuning, K e is the total taper-cavity coupling rate, K = K, ■ + K e is the total cavity decay 
rate, with K - , the intrinsic cavity damping rate, and aj n is the taper input field, which together with the output field a Q ut obeys the 
boundary condition 



^in 



(A8) 



The last two terms on the right-hand-side of eq. ( A7 1 represent the vacuum inputs due to coupling with the intrinsic (loss) bath 



of the cavity and the backward fiber taper waveguide mode, respectively (these input terms are ignored going forward as they are 



da , 



in the vacuum state and do not modify the classical field equations). In steady state, where 3? = 0, the intracavity field operator 



fl() : 



2 JA+f ' 



flin is normalized to the power incident on the cavity P ln as P; n = h(£)\(a[ n a m ) such that the intracavity photon number is 
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Combining eq. ( A8 i with eq. ( A9 1 yields the intensity transmission function 



T(A) 
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shows an example transmission 



This function describes a Lorentzian absorption curve that dips to T c { = -\ at A = 0. Figure 

curve of the device studied in this work obtained by scanning an external cavity diode laser across the fundamental resonance of 

the zipper cavity. The slope of the curve is given by 
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Usually, we lock the probe laser to a red-side detuning of A = — Jc/2, where the transduction is maximum for fixed « cav . At that 
detuning, the intracavity photon number is given by 



w cav,K:/2 



= 1- 
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(A13) 



and the slope of the transmission curve is 
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Appendix E: Derivation of the optomechanical accelerometer transduction 



Our device operates deep in the sideband unresolved regime, where co m <C K (co m = 2% x 27.5 kHz, K = (Q c /Q = 129 GHz). 
In this regime, the intra-cavity field and hence the field transmitted through the cavity adiabatically follow changes in laser-cavity 
detuning A = ©/ — CO c created by mechanical motion of the test-mass, A = goMX- In order to calculate the optical transmission 
change Ar induced by a shift of the cavity resonance frequency A, we can therefore approximate 



AT 



AT 

aK l 



such that the frequency component of the transmitted optical power arising from a displacement x(co) is given by 
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(A16) 



(A17) 



where P m is the input power in the fiber taper waveguide at the zipper cavity and Tj; n quantifies the optical loss in the fiber taper 
waveguide between the cavity and the detector via f] m = Pdet/fln, where Pdet is the optical power reaching the detector. This 
formula relates frequency components of the transmitted optical power modulation to the mechanical motion of the test-mass. 



With eq. ( A13 1 and eq. ( A15 i, this becomes 



P m {A=K/2) = {l 



T d )—gOMllmPmX- 
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(A18) 



This optical power is measured on a Newport 2117 balanced photo-detector with switchable transimpedance gain (in these 
experiments we use g t \ = 49,600 V/W), generating a voltage output of V m ~ gtiP m - An electronic spectrum analyzer (ESA) 
calculates the electrical power spectral density of this optical sideband in units of V^/Z with Z = 50 Q. and expresses it in 
dBm/Hz. The conversion follows the relation 



PSD ESA M = 10-log 



(g«Pm(CD)Y 



•1,000 



(A19) 



Careful calibration of the parameters in eq. ( |A18 i and eq. ( A19 1 as well as the optical input power, allows one to calculate the 
optomechanical coupling goM from the magnitude of the (known) thermal Brownian motion noise of the mechanical oscillator. 
In the measurements presented in Figs.[2]and[3] we have T ( / = 0.87, Q = 9,500, C0 o — 2nx 195 THz, and Tj; n = 0.57. At low 
optical input power, where negligible back-action cooling is being performed on the fundamental in-plane mechanical mode of 
the suspended test mass and the mode's effective temperature is the temperature of the room temperature bath (T ~ 3QQK), the 
optomechanical coupling constant is estimated to be goM = 2ft x 5.5 GHz/nm from the area under the Lorentzian centered at 
27 kHz of the optically transduced displacement noise PSD. This corresponds to an optical displacement sensitivity of P m /x — 
3.7 nW/pm for the fundamental in-plane mechanical mode of the suspended test mass. From electromagnetic finite-elements 
simulations we calculate goM = 2nx 13.5 GHz/nm for dimensions of the zipper cavity as measured with a scanning electron 
microscope, in good agreement with the measured value. 



Appendix F: Acceleration sensitivity measurement 



For applying AC accelerations to our device, we constructed a shake table comprising a sample holder plate glued on a shear 
piezo actuator. Applying a sinusoidal AC-voltage to the piezo creates a displacement XQsin(cOdt), which results in an applied 
acceleration of — xq(Oj sin((O c it) . For calibration of the shake table assembly, we use commercial accelerometers from Analog 
Devices of 5.5 kHz (ADXL103) and 22 kHz (ADXL001) bandwidth, respectively. In order to measure the frequency response 
of our optomechanical accelerometer, we apply a constant-voltage drive to the piezo and tune its frequency, while measuring the 
photodetector output on a lock-in amplifier. After normalizing for the o ( / -dependence of the applied acceleration, this yields the 
frequency-dependent sensitivity of the device. Normalizing an optical noise PSD then allows us to calibrate the noise-floor of 
the accelerometer in terms of a noise-equivalent acceleration. 
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FIG. A2: Autocorrelation trace of the thermal noise driven mechanical amplitude. The signal was obtained from computing the auto- 
correlation of the slowly varying magnitude of the mechanical motion returned from a lock-in amplifier. Fitting the trace with an exponential 
decay yields the time constant and thereby the mechanical (9-factor of the mode (Q m = 1 .4 x 10 ). 



Appendix G: Mechanical spectroscopy and autocorrelation method to determine mechanical quality factor 



Motion of the mechanical oscillator results in amplitude-modulation of laser light transmitted through the fiber taper which 
can be measured by monitoring the power spectral density of the detected balanced-photodiode photocurrent on an electronic 
spectrum analyzer (ESA) from which we can extract the resonance frequency and the total power in transduced sideband (pro- 
portional to mode temperature). However, the sub-Hz linewidths of our mechanical modes make establishing the quality factor 
from a measurement of the power spectral density on a spectrum analyzer infeasible because it requires a fractional stability of 
the frequency to greater than S> l/Q m ~ 5 x 10~ 7 over a period much longer than the decay time Q m /(Q m sw 10 s. To overcome 
this limitation, we extract Q„, from the autocorrelation function of the mechanical motion |30|. Since the system is driven by a 
Gaussian thermal noise process, the autocorrelation of the amplitude (X(t)X(t + %)) can be shown to decay as e~ r ' T from which 
the quality factor can be obtained as Q m = tco,„ |30|. The slowly-varying envelope of (X(t)) is obtained from the magnitude 
channel of a lock-in amplifier tuned to the mechanical resonance frequency with a bandwidth (f=a 100 Hz) much larger than the 
linewidth which ensures that small frequency diffusion does not affect the measurement of the envelope. To obtain the bare 
mechanical Q-factors the measurement is made at an optical power low enough to ensure there is no backaction. The autocor- 



relation is numerically computed and the decay is fit to an exponential curve with a constant (noise) offset. In Fig. A2 we show 
an autocorrelation trace of the device studied in this work calculated from s=s 3000 s of data sampled at 100 Hz and fit it to find 
T = 7.85 s and for (O m = 2n x 27.5 kHz that Q = 1.4 x 10 6 . For lower-Q structures, it was confirmed that this technique agrees 
with a direct measurement of the linewidth from a spectrum analyzer. In order to avoid air-damping, measurements are carried 
out in vacuum. 



Appendix H: Optomechanical and thermo-optical backaction 



The relatively small test mass makes the device studied in this work highly susceptible to optomechanical and thermo-optical 
back-action effects. Such dispersive couplings are well known to renormalize the frequency and damping rate of the mechanical 
oscillator. In particular, thermo-optical coupling which arises from a refractive index change of the material upon the absorption 
of cavity photons plays a significant role in these devices due to the efficient thermal isolation of our nano-tethered test-masses 
in vacuum. Previous studies have shown strong modification of the optomechanical spring effect and damping in similar zipper 
cavity devices ifTUll . 

The Supplementary Information of Ref. flOl gives a detailed derivation of the renormalized oscillator frequency and damping 
rate under the influence of optomechanical and thermo-optical coupling. The system of differential equations that describes the 
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time evolution of the intra-cavity field a, the oscillator position x, and the cavity temperature shift Ar is given by 



a = -[iA-(g MX + g 6l AT)]a--a- 

„• ,J2 %OM | ,2 

x = -yx-co m x \a\ 



AT = -y th AT + K abs c th hco c \a\ 



(A20) 

(A21) 
(A22) 



where g t h = — (dn/dT){co c /n) is the thermo-optical tuning coefficient, dn/dT is the thermo-optic coefficient of the material, K a ], s 
is the optical loss rate due to material absorption, c t h is the thermal heat capacity, and Yth is tne decay rate of the temperature. 
Linearizing these equations yields the static solutions 
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(A23) 



with the renormalized detuning A' = A — goMXQ — gthATb arising from the static optomechanical and thermo-optical shift. Using 
a perturbation ansatz x(t) = xq + ecos(<B m f) one arrives after some algebraic manipulation at a modified harmonic oscillator 
equation for x with a renormalized frequency co' m and damping rate / given by 
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where the transfer function g(co) is defined as 
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and Ath = gthATb is the static thermo-optical shift of the cavity resonance frequency. In the sideband unresolved regime where 
CO m <^.K and for thermal decay rates Yth smaller than the mechanical frequency, an approximation of g(a>) yields 
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with the correction factors 
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and the saturation parameter 
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FIG. A3: Demonstration of thermo-optomechanical damping and cooling. The green bullets show measured g-factors of the mechanical 
mode as function of the optical power, yielding thermo-optomechanical damping by a factor of ~ 280. The blue bullets show the corresponding 
optical power in the sideband generated by mechanical motion, proportional to the effective mode temperature. We observe cooling to 
T e g ~1K. The dashed green curve corresponds to a theoretical model that includes optomechanical and thermo-optical back-action. 



In the parameter regime of our devices, purely optomechanical back-action is a relatively weak effect due to the low optical Q- 
factor. For the parameters given above and for a pump laser with an incident power of P m = 116 /iW half a linewidth red-detuned 
from the cavity resonance, optomechanical back-action alone predicts a frequency shift of merely co' m — CO,„ = — 2n x 35.9 Hz 
and a damping factor of Y m /y = 1.01. 

In order to study the influence of thermo-optical back-action, we measured the Q-factor of the mechanical mode as function 
of the optical power launched into the cavity, shown as the green bullets in Fig. A3 When increasing the optical power to 
Pi n sa 300 jtiW, which corresponds to an intracavity photon number of n cav ~ 1, 100, the g-factor shows strong damping and is 
reduced by a factor of w 200. Similarly, we measure the area of the mechanical resonance peak from the optically transduced 
thermal noise PSD for a series of optical powers, and plot the inferred effective mode temperature as blue bullets in Fig. |A3| 



Clear in Fig. A3 is that the effective mode temperature is dropping with the measured mechanical g-factor. 



The observed mechanical damping is much larger than the value predicted by pure optomechanical back-action and can be 



explained when includi ng the rmo-optical tuning. The green line in Fig. A3 was obtained by calculating the modified g-factor 
Q' m = (Om/Y,„ using eq. (A25 i with A tn = — 0.05fC and y tn = 2n x 9.2 kHz. The latter value is in good agreement with the one 
from flOl (7 t h = 2n x 10 kHz), which suggests that the time constant of thermo-optical tuning is dominated by heat-flow from 
the zipper cavity region to the reservoir formed by the test-mass (or the bulk in the case of flOl . respectively). 

The obtained values for A tn and y tn result in correction factors of V = 12,400, W = —0.011, and a saturation parameter 
of s « 3 x 10~ 36 . Accordingly, we expect a significant thermo-optical correction to damping, as observed, but only a minor 
modification of the optomechanical spring: co' — CO = 2n x 36.2 Hz for the pump power used in the experiment. Indeed, we 
observed a frequency shift of 101 Hz, in reasonable agreement with the theoretical value. 



Appendix I: Analysis of optical noise power spectral densities 



As discussed above, noise power-spectral-densities (PSDs), such as those shown in Fig. [2p, arise from the contributions of 
various noise sources. In the following we derive expressions for their magnitudes. Throughout the analysis below we work 
with single-sided PSDs, unless otherwise stated, as these are the PSDs measured in our experiment. 
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1. Noise from thermal Brownian motion 

In contact with a heat-bath at room temperature, the test-mass oscillator is subjected to thermal Brownian motion. From the 
equipartition theorem, the root-mean-square displacement of a harmonic oscillator is given by 



*£. (A35) 

k 



If we assume the acceleration-noise exerted by the bath to be white, i.e. frequency-independent, its power-spectral density S' c , a 
has to obey 

*L= flxWl 2 St* da, (A36) 

Jo 

such that thermal test-mass motion corresponds to a noise-equivalent acceleration (NEA) of 



f— 4k B Tco m 4k B Tr„, 

«th = V 5 ™ = \ n = V ■ (A37) 
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In the device presented in this work, we have (0 m = 2k x 27.5 kHz , m = 10~ n kg, Q m = 1.4 x 10 6 , T = 295 K, and therefore 
flth = 1.4 /xg/y/Hz. For a mass-on-a-spring oscillator with »,„ = \Jkjm this corresponds to 

/ 4 ^^' /4 , A1 « 

y Qm ml* 

Driving the harmonic oscillator with susceptibility #(©), this NEA translates into frequency-dependent displacement noise 
according to 

S f xi( o) = 4 J^ l - 2 . (A39) 
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According to eq. ( A17 I, the optical signal transduced by the cavity then exhibits a noise power-spectral density of 
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Under the influence of thermo-optomechanical back-action discussed in appendixfH] the dynamic parameters co m and Q m have 
to be replaced by the renormalized values co' m and Q' m . With the parameter values realized in this experiment, the optical noise 

arising from thermal Brownian motion corresponds to yS'p P (Q) = 0.96 pW/-\/rTz. 

2. Optical shot noise 

Photon shot noise arises from the quantum nature of light and from the destructive character of optical measurements using 
photodiodes. The single-sided shot-noise power-spectral-density for light of frequency C0q and power Pd et incident on a photo- 
detector is frequency-independent and given by 

4?=^, (A42) 

*7qe 

where the quantum efficiency T7 qe (=0.84) is linked to the photodiode responsivity R (=1 A/W) via 

R=P*. (A43) 

h(0 o 
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In our balanced detection scheme, we consider the shot noise of the difference photocurrent of the two detectors. Since photon 
annihilation at the two detectors is uncorrelated, the total shot noise is given by the incoherent sum of the two individual power- 
spectral-densities, such that 

S™ = 2hC0 ^ (A44) 

»7qe 

with P t ot = ^deti +Pdet2 being the sum of the individual powers hitting the two photodiodes. In our balanced detection scheme, 
^deti = Pdet2 an d Ptot = 27'tJjii/V While the balanced detection scheme used in our experiment is beneficial towards the suppres- 
sion of technical laser amplitude noise, it hence comes with the disadvantage of introducing more shot noise into the system. 
In this experiment, the noise-equivalent power corresponding to shot noise is 6.1 pW/\/Hz. The noise-equivalent acceleration 
corresponding to this noise background is given by 

a SN (co) = Js™ = l -j— JS™ (A4S) 
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With the values given above, this yields «sn = 8.9 jUg/\/Hz around DC. 

3. Detector noise 

The electronic detector noise is usually quantified by the noise-equivalent-power (NEP), which for the Newport 2117 detector 
and the transimpedance gain setting we use is on the order of 2.8 pW/VHz. The optical noise power-spectral-density then is 

Spp P (co) =NEP 2 . (A47) 

In analogy to eq. ( |A46| >, the NEA corresponding to electronic detector noise can be derived as 

a N EP = 7 r— i , m NEP. (A48) 
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4. Backaction noise 



Here, this is found to be anep =4.1 jug/ yHz at DC 



The extra noise a a( jd added by the optical field mentioned above arises from optical noise that exerts a random force on the 
mechanical oscillator via radiation pressure. The optical noise arises from classical amplitude noise and from intrinsic shot 
noise. In the following, we consider only quantum back-action noise aba arising optical shot noise. With ftgoM being the force 
exerted per photon and for n cav photons in the cavity, the random acceleration created by optomechanical back-action has a 
power spectral density of [ 22 1 



aa — ^ 2 Wcav -- ' (A4y) 






resulting in a noise-equivalent acceleration of aba = V*^? = 5-6 ng/v^Hz. Here, owing to the low quality factor of the optical 
cavity and the low mechanical frequency, the shot noise radiation pressure force is approximately white noise for frequencies 
of relevance near the mechanical frequency. Note also that we are using single-sided PSDs, hence double the value of the 
(approximately) symmetric double-sided PSD. This value is much smaller than the acceleration noise created by the other 
sources discussed previously. The frequency-dependent displacement noise created by quantum back-action is 

cBA,„\ _<, (2hgOM\ "cav, ,,„\<2 
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On the mechanical resonance, and using eq. ( A13 1, this yields 
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FIG. A4: Frequency response of acceleration noise and power-scaling, a, Frequency dependent contributions of different noise sources to 
the noise-equivalent-acceleration (NEA). For this calculation, we used the same parameters as for the device shown in Fig. [2] b, Dependence 
of the DC-NEAs as function of incident laser power. The dashed vertical line indicates the optical power used in measurements of Figs.pjs 
and[3] 



resulting in y/S^((O m ) = 23 fm/\/Hz for the device and experimental conditions described in Figs.pl and p|(P m = 116 /xW). 
This should be compared to the fundamental standard quantum limited displacement noise given by 
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which on resonance has the simple form 



with the zero-point motion given by 
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For the device and experimental conditions described in Figs. 2 and 3 this yields an on-resonance SQL of ySjc =2.8 fm/\/Hz. 



5. General discussion 



The dashed lines in Fig. [2b show the contributions of these noise terms to the PSD of the balanced photo-detector output, 



where we neglected back-action noise. Figure A4 1 shows the corresponding frequency-dependent noise-equivalent acceleration 
values corresponding to the different noise terms for the device studied here. Here, we include oba, which can be seen to only 
contribute negligibly to the NEA of the device. While a± (green) and aba (gold) are frequency-independent, the NEAs of photon 
shot noise «sn (red) and electronic detector noise «nep (cyan) are colored by the frequency-dependent response of the oscillator 

%{(0). 

While thermal noise arises as a fundamental property of a mechanical oscillator in contact with a heat bath at temperature T, 
the contri bution s of shot noise and detector noise are dependent on the efficien cy of t he optomechanical transduction mechanism. 



From eq. ( A46 1 and eq. ( A44 1, one can see that osn x P\ 



-1/2 



,1/2 



, while from eq. (|A48|) it follows that anep x P^ ■ Similarly, back 

For illustration, Fig. 



A4 3 shows the 



action noise scales with the square-root of the number of photons in the cavity: aba x P\ 
relative contributions of the individual noise terms at DC (CO = 0) as function of incident power. Here, we include the effects of 
thermo-optomechanical back-action on the mechanical susceptibility, as discu ssed in appendixlH] The thermal noise background 
is not affected by cooling of the mechanical mode, since it follows from eq. ( A37 i that aj? h °< r e ff y m — const under back-action 
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damping/cooling of the mechanical mode. The roll-off of shot noise and detector noise for pump powers above 10 mW arises 
from the decrease of the mechanical mode frequency co'„ due to the optomechanical spring effect. This results in an increase of 
the DC acceleration sensitivity j(0) = l/<0,'„ and thereby a reduction of the corresponding acceleration noise floors according 
to eqs. ( |A46| ) and ( |A48| >. 

As mentioned previously, for the power used in the experiment, the NEA is limited by photon shot noise. For two orders 
of magnitude higher pump powers, the NEA starts being dominated by thermal noise of the test-mass oscillator. Alternatively, 



according to eq. ( A46 1, thermal-noise limited detection can be achieved by increasing goM by one order of magnitude. 



Appendix J: Linear Dynamic Range 



A key requirement for any inertial sensor is linear response over a reasonable dynamic range. To check the linearity of the 
response of the accelerometer presented in the text, we varied the amplitude of a sinusoidal signal sent to the shear piezo at 
9.92 kHz and recorded the voltage corresponding to the peak height of the transduced modulation tone - shown in blue bullets 
in Fig. A5 The sensor behaves linearly over a dynamic range of 41 dB, with the tone vanishing into the shot noise floor for 



an applied acceleration of « 10 jig at a resolution bandwidth of 1 Hz. The green bullets in Fig. A5 show data from a different 
device with similar geometry but slightly lower mechanical Q-factor, which exhibits a linear response over 49 dB. This particular 
measurement was limited by the maximum output voltage of the function generator. Ultimately, however, the linear dynamic 
range ends when motion of the test mass shifts the optical resonance by a magnitude comparable to the optical cavity linewidth. 
For this device, this is expected to occur for accelerations of ~ 50 g for frequencies below ~ 25 kHz. 
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FIG. A5 : Demonstration of large linear dynamic range of typical devices. While varying the amplitude of the acceleration applied with the 
calibrated shake table at 9.92 kHz, we measure the optical signal transduced via the mechanical mode. The blue bullets show the transduced 
signal of the device presented in the text using the voltage corresponding to the peak height of the modulation tone on the ESA spectrum, 
which exhibits linear response over 40 dB. The inset shows the corresponding PSD spectra from the ESA for modulation tones between 0.1 g 
and 12.8 fig, taken at a resolution bandwidth of 1 Hz. The green bullets show data obtained from a different device with a larger thermal noise 
background but very similar optomechanical coupling using the lock-in scheme depicted in Fig.pt, which exhibits linear response over 49 dB 
- limited by the maximum voltage output of our function generator, which corresponds to an acceleration of 8 g. The black lines are linear fits 
to the data. 



